Split of zero-bias conductance peak in normal-metal / d wave superconductor 

junctions 
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Effects of impurity scatterings on tlie conductance in normai-metal / d wave superconductor 
junctions are discussed by using tlie singie-site approximation. So far, tlie split of tlie zero-bias 
conductance peak lias been believed to be an evidence of the broken time reversal symmetry states 
at the surface of high-Tc superconductors. In this paper, however, it is shown that the impurity 
scattering near the interface also causes the split of the zero-bias conductance peak. Typical con- 
ductance spectra observed in experiments at finite temperatures and under external magnetic fields 
are explained well by the present theory. 

PACS numbers: 74.81.-g, 74.25. Fy, 74.50.-|-r 



I. INTRODUCTION 

The zero-energy state (ZES)! formed at surfaces of 
superconductors is a consequence of the unconventional 
symmetry of Cooper pairs. Since the ZES appears 
just on the Fermi energy, it drastically affects trans- 
port properties through the interface of junctions con- 
sist of unconventional superconductors For instance in 
normal-metal / high- Tc superconductor junctions, a large 
peak is observed in the differential conductance at the 
zero bias voltagei '^'^i^i^i''i^i^i^° The ZES is also responsi- 
ble for the low-temperature anomaly of the Josephson 
current between the two unconventional superconduc- 
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An electron incident into a normal-metal / supercon- 
ductor (NS) interface suffers the Andreev reflection^^ by 
the pair potential in the superconductor. As a result, a 
hole traces back the original propagation path of the inci- 
dent electron. This is called the retro property of a quasi- 
particle which supports the formation of the ZES. Strictly 
speaking, the electron-hole pairs just on the Fermi energy 
hold the retro property in the presence of the time rever- 
sal symmetry (TRS). Thus the ZES is sensitive to the 
TRS of the system. Actually, the zero-bias conductance 
peak (ZBCP) in NS junctions splits into two peaks un- 
der magnetic Helds.^^'^^t^^^ The peak splitting is also 
discussedSLSSiSSi^fliSiiSiS^ when the broken time reversal 
symmetry state (BTRSS) is formed at the interface. The- 
oretical studies showed that such BTRSS's are character- 
ized by the s-\-idx^^ or dxy-\-idx2_y7^ wave pairing sym- 
metry. Experimental results, however, are still controver- 
sial. Some experiments reported the split of the ZBCP at 
the zero magnetic field)'^^i'^^i^'^i'^^i^^i^°i^^i'^^ other did not 
observe the splitting.SiS42i4^ The ZBCP is also sensi- 
tive to the exchange potential in ferromagnets attaching 



to unconventional superconductors 4^*^ 

In previous papers, we numerically showed that ran- 
dom potentials at the NS interface cause the split of the 
ZBCP at the zero magnetic field by using the recursive 
Green function method. We also showed that the 
splitting due to the impurity scattering can be seen more 
clearly when realistic electronic structures of high- Tc ma- 
terials are taken into account?^ Unfortunately, we could 
not make clear a mechanism of splitting. Our conclu- 
sion, however, contradicts to those of a number of theo- 
ries^ii^SiSi^^iSiS based on the quasiclassical Green func- 
tion methodiSSiSSiSiSSiSi The drastic suppression of the 
ZBCP by the interfacial randomness is the common con- 
clusion of all the theories. The theories of the quasiclassi- 
cal Green function method, however, concluded that the 
random potentials do not split the ZBCP. Thus this issue 
has not been fixed yet. There are mainly two reasons for 
the disagreement in the two theoretical approaches (i.e., 
the recursive Green function method and the quasiclassi- 
cal Green function method) . One is the treatment of the 
random potentials, the other is the effects of the rapidly 
oscillating wave functions on the conductance. In our 
simulations, we calculate the conductance without any 
approximation to the random potentials and the wave 
functions; this is an advantage of the recursive Green 
function method42iS^ 

In this paper, we discuss effects of the impurity scat- 
tering on the conductance in normal-metal/d wave super- 
conductor junctions by using the Lippmann-Schwinger 
equation. We assume that impurities are near the NS in- 
terface on the superconductor side. The differential con- 
ductance is analytically calculated within the single-site 
approximation based on the conductance formula^^*^ 
The split of the ZBCP due to the impurity scattering is 
the main conclusion of this paper. The impurity scatter- 
ing affects the conductance in two ways: (i) drastically 
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suppressing the conductance around the zero bias voltage 
and (ii) making the conductance peak wider. The spht 
of the ZBCP is a consequence of the interplay between 
the two effects. In the present theory, we successfully 
explain typical conductance shapes observed in several 
experiments. We also show that the splitting peaks are 
merged into a single conductance peak for sufficiently 
high temperatures and that the peak splitting width in- 
creases with increasing external magnetic fields. 

This paper is organized as follows. In Sec. II, we 
derive the reflection coefficients in NS junctions within 
the single-site approximation based on the Lippmann- 
Schwinger equation. The split of the ZBCP is discussed 
in Sec. III. In Sec. IV, calculated results are compared 
with experiments and another theories. In Sec. V, we 
summarize this paper. 



II. LIPPMANN-SCHWINGER EQUATION 

Let us consider two-dimensional NS junctions as shown 
in Fig.^ where normal metals (x < 0) and d wave super- 
conductors {x > 0) are separated by the potential barrier 
VB{r) = Vo6{x). We assume the periodic boundary con- 
dition in the y direction and the width of the junction 
is W. The a axis of high-Tc superconductors is oriented 
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FIG. 1: The normal-metal / d-wave superconductor junction 
is schematically illustrated. The crosses represent impurities. 

by the 45 degrees from the interface normal. The pair 
potential of a high- Tc superconductor is described by 



(1) 



in the momentum space, where Aq is the amplitude of 
the pair potential at the zero temperature, = cos 7 = 
kx/kp and ky = sin 7 = ky/kp are the normalized wave 
number on the Fermi surface in the x and the y direc- 
tion, respectively. The Fermi wave number kp satisfies 
h^kp/2m = fip, where fip is the Fermi energy. The 
schematic figure of the pair potential is shown in Fig. ^ 
The NS junctions are described by the Bogoliubov-de 



Gennes equation^SSi 
dr' 
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where Rc — {Xc,Yc) = (r + r')/2 and = r — r'. 
Throughout this paper, we neglect the spatial depen- 
dence of the pair potential near the junction interface. 
This is a reasonable approximation when we consider the 
conductance around the zero bias voltage?^ The spatial 
dependence of the pair potential should be determined in 
a self-consistent way when we discussed the conductance 
far from the zero bias such as eV ~ Aq. Here V is the 
bias voltage applied to junctions. 

We consider impurities near the interface on the su- 
perconductor side as indicated by crosses in Fig. ^ The 
potential of impurities is given by 



Viir) = V,J2S{r-r,), 



(6) 



where Ni is the number of impurities. In the absence 
of impurities, the transmission and the reflection coeffi- 
cients are calculated from boundary conditions of wave 
functions at the junction interface as shown in Ap- 
pendix A. By using these coefflcients, four retarded Green 
functions are obtained as shown in Appendix B. The nor- 
mal conductance of the junction is given by 



h 

cos^7 
Jo Zq + cos^ 7 



(7) 

(8) 



where Tb is the transmission probability of the junction, 
Nc — 2W/Xp is the number of the propagating channels 
on the Fermi surface, Xp — 2'iT/kp is the Fermi wave 
length and zq = mVo/{h^kp) represents the strength of 
the potential barrier at the NS interface. In the limit of 
Zq >• 1, Tb is proportional to l/zg. 

Effects of impurities on the wave functions are taken 
into account by using the Lippmann-Schwinger equation, 

^'^'\r) =4'\r) + jdr'G^{ry) Vi{r') a, iji'\r'), (9) 

Ni 

=i^l'\r) + J2Goir,rj) V, ^^^rj), (10) 
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where I indicates a propagating channel characterized by 



the transverse wave number k 



Here ipg' 



(r) is the wave 
(0 



function in which an electronhke quasiparticle with ky is 
incident into the clean NS interface from normal metals 
and is described as 
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for a; < 0, where = \/kf ± kpE / jip is the wave num- 
ber of a quasiparticle in normal metals and k^ + fcy ^ = 
k'jp. For a; > 0, the wave function in clean junctions is 
given by 
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where is a macroscopic phase of a superconductor, 
k^ — (kf ± kj^fli/ fip^^^^ is the wave number of a quasi- 
particle in superconductors, il; = y/ E'^ — Af and A; — 

2AQkiky\ The wave function at an impurity ijA''\rji) 
can be obtaind by r ^ r^v in Eq. (fTn|l 



1 



^(')(r,). 



(17) 



It is possible to calculate the exact conductance if we ob- 
tain ip^^^rj) for all impurities by solving Eq. p7|) . Ac- 
tually it was confirmed that the conductance calculated 
from the numerical solution of Eq. (|17|l is exactly identi- 
cal to that computed in another numerical methods such 
as the recursive Green function methodnSiSl In this pa- 
per, we solve Eq. H17|l within the single-site approxima- 
tion, where the multiple scattering effect involving many 
impurities (Anderson localization) are neglected. How- 
ever the multiple scattering by an impurity is taken into 
account up to the infinite order of the scattering events. 



In the summation of j in Eq. (|17|l . only the contribu- 
tion with j = j' is taken into account in the single-site 
approximation^ In this way, the wave function at rj is 
approximately given by 



SSi 



4'\r,). (18) 



(To -G^''(rj,rj) Vi aa 

We note that the single-site approximation yields the ex- 
act conductance when Ni = 1. Within the single-site 
approximation, Eq. ffTU\i can be solved as 



rsUr) =4" (r) + J2 ^o'(.r, r,) V, a. 



X <7o-Gf'(^j-,rj)^.^3j (19) 

for a; < 0. On the right hand side of Eq. (|19|) . all functions 
have been given by analytical expressions. 

In the presence of the impurity scattering, the wave 
function Eq. (fT^ can be expressed as 



in . X A I 

e Ai' I + 
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for X < 0, where Aii^i and Bii^i are the Andreev and the 
normal reflection coefficients in the presence of impuri- 
ties, respectively. These coefficients are obtained from 
relations 
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The scattering theory based on the Lippmann-Schwinger 
equation requires complicated algebra as shown below 
because the perturbation expansion is carried out in the 
real space. In return, effects of the impurity scattering 
can be taken into account up to the infinite order of the 
perturbation expansion without using any self-consistent 
treatments. In addition, the reflection coefficients are 
explicit functions of the impurity positions in a single 
disordered sample. These are advantages of the present 
method. 

In what follows, we consider low transparent junctions, 



(i.e. 



^ 1). From the reflection coefficients in Ap- 



pendix A, the Green function in the superconductor is 
given by 
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where p : 
defined by 



n 



2 cos 7 /ii; 



The local density of states^ at r is 



Ns{E,x) 



--ImTrG^'^(r,r) 



(24) 



The first term in Eq. (|23|l contributes to the bulk den- 
sity of states. Since 2p is roughly estimated to be i/£,o, 
another terms contribute to the local density of states 
near the interface, where ^0 — ^ivf/ttAo is the coherence 
length and vp = hkp/m is the Fermi velocity. In low 
transparent junctions, 6th, 7th and 8th terms are neg- 
ligible. The 4th and the 5th terms are also negligible 
because integrals of such rapidly oscillating functions be- 
come very small. The 2nd and 3rd terms are dominant for 
E ^ Aq. In Fig. 121 (a), we show the trace of the Green 
function in Eq. (|23() as a function of E, where 20=10, 
Ao = 0.1/ii?, xkp = 6, and Nq = m/(7r?i^) is the normal 
density of states in the unit area. Since ^o^f ~ 6.3, the 
results correspond to the Green function at a distance ^0 
away from the interface. The horizontal axis in Fig.|2(a) 
is normalized by 



EzEP — 



Ao 




(25) 
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FIG. 2: In (a), the trace of the Green function in the su- 
perconductor is shown as a function of E, where 20 ~ 10, 
xkp = 6 and Ao = O.lfip- The peak width of the imaginary 
part is given by Ezep = Aq/zq. The energy scale, -Edip, char- 
acterizes the drastic increase of the imaginary part and the 
drastic decrease of the real part. In (b), the local density of 
states is shown as a function of xkp, where E = 0, zo = 10 
and Ao = O.I/if. The numerical and the analytical results 
are denoted by the solid and the broken lines, respectively. 



The solid and the broken lines represent negative of the 
imaginary part and the real part of Eq. (|23|) , respectively. 
As shown in Fig. [3 (a), the imaginary part of the Green 
function has a large peak around E — reflecting the 
ZES formed at the junction interface. The energy scale, 
Ezep, characterizes the width of the zero-energy peak. 
The real part of the Green function first increases with 
decreasing E then suddenly decreases to zero for £' ^ 0. 
The detail analysis indicates that the real part of the 
Green function has its maximum around an energy 



E, 



dip 



Ao 



Zo(xfcF)'^ 



(26) 



The Green function for E < Ezep is approximately cal- 
culated from the 2nd and the 3rd terms of Eq. (|23() 

G^^(r,r) « 27rNoz'oe--/i°ig2 ~ tgi)ao, (27) 
2 f"^^ Aq cos*7 sin^7 sin^(a;A:F COS7) 



91 



c?7 



(28) 



TTJo E^ + Al cos^ -y s'm^ ^ 

2 /"^^j £^Zq Ao C0S7 sin7 sin^(a;fcF COS7) 
TT Jo E^Zq + Aq cos^ 7 sm 7 



where we use fl ~ i|2AoCOS7sin7|. The imaginary part 
of the Green function, gi , is of the order of unity when 
E ~ Ezep- However, gi at E = becomes much larger 
than unity for xkp ^ ^ because 



2 sin2{x/cFCOS7} 
gi{E = 0) = - / d7 ■ 

71" Jo 



cos^ 7 



xkf 



(30) 



Thus the energy scale, E'dip, characterizes the drastic in- 
crease of gi and the drastic decrease of g2- The local 
density of states at £^ = calculated from Eqs. H23|l and 
(|24|1 is plotted as a function of xkp in Fig.[2Ib) with the 
solid line. For comparison, we also show the analytical 
results represented by 



2TrN( 



E=0 



XkF<Jo, 



(31) 



with the broken line. The results show the remarkable 
enhancement of the local density of states around a; ^ ^o- 
This implies that the ZES is formed around x ^ S^o- 

Here we note several remarks as follows. To calculate 
the Green function, we consider the 3rd term in Ea. H23|) 
which rapidly oscillates as cos(2a;fcF C0S7). Such rapidly 
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oscillating terms are usually neglected in the quasiclassi- 
cal Green function method. We, however, cannot neglect 
the 3rd term because it removes the divergence of the 
local density of state dX E = Oi^Siii The 3rd term also 
becomes important when we calculate the local density 
of states just at the surface, (i.e., x — 0), 
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No TT Jo 
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2E^zl 
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2 /Ao 

-K \ E 



2 r^^i^^zo +^ocos^7sin^7 
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(32) 
(33) 



where K{x) is the complete elliptic integral of the first 
kind and describes the bulk density of states. Another 
terms come from the 4th and 6th terms in Eq. H23I) . The 
first equation is the exact expression and we use E < 
EzEP in the second line. We exactly obtain N{E = Q,x ^ 
0) = A^o- Thus there is no remarkable enhancement in the 
zero energy local density of states just at the interface. 
The 2nd and the 3rd terms in Eq. do not contribute 
to the Ns{E,0) because they exactly cancel with each 
other at a; = 0. 

In the next section, the conductance for E < E^ep will 
be discussed. It is possible to rewrite a part of Eq. H19|) 
as 



si = Qo gi, 

S2 = Qo 92, 
S = S2~ isi, 



l-s2 ' 



(34) 

(35) 
(36) 
(37) 
(38) 



where s corresponds to the self-energy of a quasiparticle 
scattered by an impurity once and describes two impor- 
tant features. First, s becomes large in low transparent 
junctions even if ViNo is fixed at a small constant because 
Qo in Eq.(|2Z|) is proportional to V^Nqz^ ~ ViNq/Tb- 
Thus Qo represents the normalized strength of the im- 
purity scattering. This behavior explains the previous 
numerical simulation.'*^ Second, effects of impurity scat- 
tering far away from the interface on the conductance is 
negligible because Qo decreases exponentially with the 
increase of x. Impurities around the ZES, (i.e., x ~ ^o) 
seriously affect the conductance for E < Ezep- 



III. CONDUCTANCE 



The differential conductance in NS junctions is calcu- 
lated from the normal and the Andreev reflection coefli- 



l,m 
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dfME-eV) 



d{eV) 

X [&un-\B„a{E)? + \An,l{E)\''] , 




(39) 
(40) 

(41) 



where fpuiE) is the Fermi-Dirac distribution function. 
When Zg ^ 1 and E < Ezep, the reflection coefficients 
are calculated as 

+e-'''^L„,i [A„|A,|(s + 1) + A,|A„,|(s - 1)] , 

(42) 

+ iL„,i[|A„||Ai|(s-(-l) + AiA™(s- 1)], (43) 

n N- 



(1 - s^)ki 



3 = 1 



^ V^ g»(p,„+P,)^. sin{k^xj) sm{kixj). (44) 

The conductance is then given by 

fdfFD{E~eV) 
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V d{eV) 
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Aq cos'^7 sin^7 



A§ COS67 sin^7 ' 



Tj =Re 
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h =- 

IT 
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sQo{2l2 - ill + ih) 







Aq COS "7 sin 7 sin {xjkp cos j) 
{E^z^ + Al cos67sin27)2 ' 

Ezq Aq cos'^7 sin^7 sin^(xj k p cos 7) 
{E^z^ + A^o cos67sin27)2 ' 
E'^ZqAI cos'*7sin^7sin^(xjfc_F C0S7) 
^ {E^z* + Al cos67sin27)2 ' 



c?7 

Tr/2 



(45) 

(46) 
(47) 
(48) 
(49) 
(50) 



The first term of Eq. ggi, A^c5^°\ is the conductance 
in clean junctions and Tj represents effects of the impu- 
rity scattering on the conductance. When we calculate 
|^m_;|2 and |_B,„_;|2, thc summation with respect to im- 
purities Sj^' must be carried out only for j' = j 
in the single-site approximationiM As a consequence, the 
current conservation law is satisfied for E < Ezep- 

To study effects of impurities on the conductance, we 
first assume Xj = xo for all impurities. The conductance 
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is rewritten as 



2e^ p ( dhT>{E-eV) 



(51) 



where rii = NiXp/W is the dimensionless hne density of 
impurities less than unity. When scattering effects are 
strong, \s\ ^ 1, Ni cannot be much larger than W/Xp- 
This limits the applicability of the single-site approxima- 
tion. 

We show conductance for several choices of xo/c^ and 
ViNfj in Fig. |3| where zq — 10 and rii — 0.9. The two 
parameters are chosen as xt^kp ~ 10, V^iVo — 0.01 in 
(a), xokp = 2.0, ViNo = 0.005 in (b), xakp = 26.0, 
V,No = 0.1 in (c) and x^kp = 12.0, V.Nq = 0.1 in (d). 
The broken line is the conductance in clean junctions. 
The temperature is fixed at a very low temperature T = 
O.Oli^ZEP which is estimated to be 0.05 K by using Aq = 
50 meV for zq = 10. As shown in (a)-(c), the ZBCP 
is splitting into two peaks by the impurity scattering. 
While the results in (d) shows the single ZBCP. Roughly 



scattering is negligible as shown for Xokp = 50 because 
Qo in Eq. (|S7|l becomes almost zero. For xgkp = 2, 5 
and 10, Tj increases with decreasing E, which indicates 
the enhancement of the impurity scattering around E — 
0. The suppression of the conductance around the zero 
bias is explained in terms of the drastic increase of the 
local density of states with decreasing E as shown in 
Fig- [21 (a)- Therefore the suppression of the zero- bias 
conductance happens irrespective of xokp and ViNo. In 
addition, a nonmonotonic E dependence of F^- for xgkp — 
10 is a source of small conductance peak at the zero-bias 
in Fig.|31(a). The same small peak is also found in Fig.|3| 
(b). 

The widening of the ZBCP can be explained by the 
sign change of Tj. When E > 0.15i?zEP in Fig. ^ F^ for 
Xokp = 10 becomes negative and impurities enhance the 
conductance. As a consequence, the conductance peak 
becomes wider than g^^\ The split of the ZBCP is a 
consequence of the interplay between the suppression of 
the conductance around the zero bias voltage and the 
widening of the ZBCP as shown in Fig.|31(a)-(c). There- 
fore the sign change of F^ explains the split of the ZBCP. 
For E(\ip < E < EzEP , Pj is almost a decreasing function 




E/E 

ZEP 



FIG. 3: The conductance is plotted as a function of bias 
voltages, where zo — 10 and rii = 0.9. The parameters are 
Xokp = 10, KNq = 0.01 in (a), xokp = 2.0, V^No = 0.005 in 
(b), xokp = 26, ViA^o = 0.1 in (c) and xokp = 12, ViNo = 0.1 
in (d). The broken lines denote the conductance in clean 
junctions. 

speaking, the impurity scattering affects the ZBCP in 
two ways: (i) it decreases the conductance around the 
zero bias voltage and (ii) it makes the ZBCP wider. The 
two effects (i) and (ii) are well characterized by the E 
dependence of Tj and the sign change of Tj in Eq. (|51|l . 
respectively. In Fig. 0] Tj is plotted as a function of 
E for several Xokp, where ViNo is fixed at 0.01. We 
note that Tj for xgkp — 10 yields the conductance in 
Fig. 121 (a). When xgkp ^q, ~ and impurity 



FIG. 4: The function Tj{E) is plotted as a function of E for 
several xokp at ViNo — 0.01. The arrows indicate iJdip. For 
xokp = 10, Edip = 0.001 EzEP. 

of E and is positive at £' = E^ip as shown in Fig. ^ It 
is followed from Eq. H47() that 

Tj oc {si(/i - h) + 2/2S2} + {.si(/i - h) - 2/2S2} . 

(52) 

Within our study, si(/i — I3) tends to be much smaller 
than 2/2 S2 for -Bdip < E < Ezep, which implies an im- 
portance of the real part of the self-energy, S2, for the 
splitting. The sign change of Tj happens when the im- 
purity scattering is sufficiently weak so that 

\s\^ = sl + sl^l, (53) 
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is satisfied. When \s{E = -Bdip)| is a small value less 
than unity, the effects of impurities are negligible and 
the conductance almost remains unchanged from that 
in clean junctions. The split also cannot be seen when 
\s{E = EzEp)\ is larger than unity. An example is 
shown in Fig. El^d), where xokp — 12.0, ViNo — 0.1 and 
\s{E = EzEp)\ is estimated to be 1.5. In this case, the 
suppression of the zero-bias conductance dominates over 
the widening of the ZBCP. As a result, the conductance 
is always smaller than that in clean junctions and the 
ZBCP remains in a single peak. 

We should pay attention to a similarity in the shapes 
of the conductance in the present theory and those in 
experiments. Amazingly, the conductance structure in 
Fig. 121 (a) is very similar to that observed in the exper- 
iment^. It is possible to find a very small conductance 
peak at 1^ = in addition to the splitting peaks around 
V ^ ±lmV in Fig. 2 of Ref. |3^ The impurities away 
from the interface explains another conductance shape 
in the experimenli^. The conductance structure in Fig. 2 
of Ref. f§ is very similar to that in Fig. |3I (d) . The present 
theory explains, at least, two typical conductance shapes 
observed in the experiments. 

As shown in Fig. |31 the magnitude of the impurity 
potential and the position of impurities are key factors 
for the degree of splitting. In Fig|Sl the gray area indi- 
cates sets of (ViNo, xokp) which satisfy Eq. H53|) within 
-Edip < E < EzEP- The open circles denote sets of 
{ViNo, xokp), where we find the split of the ZBCP in 
Eq. (|51|l . All the circles are inside of the gray area. 
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FIG. 5: A phase diagram for the split of the ZBCP. The 
gray area indicates sets of (ViA^o, xokp) which satisfy Eq. 1531 . 
The open circles denote sets of {ViNo, xokp), where we find 
the split of the ZBCP in Eq. (I51II . The local density of state 
at is = in Fig. |21 (b) is also shown. The filled squres are 
parameters used for the conductance in Fig. |3 

Although the circles and the gray region are not per- 
fectly coincide with each other, they show qualitatively 
the same tendency. The parameters used in Fig. |3| are 
indicated by filled squares. Since s at E' = is propor- 
tional to the local density of states, we also show the 



Ns{E = 0,x) in Eq. To satisfy Eq. impurities 

around xq ~ should have sufficiently small scattering 
potentials because the local density of states has large 
values there. Thus the gray area appears for small impu- 
rity potentials near the interface, (i.e., xq ^ ^o)- The gray 
region spreads to larger ViNq as the increase of xq > 
because the local density of states becomes smaller val- 
ues. The results imply that the strong impurities are not 
necessary for the split of the ZBCP. It is evident that 
this phase diagram is valid in the limit of high impurity 
density and the diagram would be changed depending on 
the transmission probability of junctions. 

In Fig. O all impurities are aligned at Xj = xq. In 
real junctions, however, impurities may be distributed 
randomly near the interface as shown in Fig. ^ The con- 
ductance in such realistic junctions are shown in Fig. El 
where impurities are distributed randomly in the range 
of 1 < Xjkp < Lgkp, Pi = Ni ^£ is the dimensionless 
area density of impurities and zq — 10. The conductance 
is calculated from an expression 



(r) 



(54) 
(55) 



i=i 



where (• • • ) represents the ensemble average. Since the 
conductance in Eq. H45|) is characterized by the number 
of impurities, a factor kpLs/ir appears in Eq. H54|l . We 
choose Lgkp = 20 in Fig. Elbecause we focus on the im- 
purities near the interface and ^q^f ~ 6.3. We consider 
low density strong impurities in (a), where ViN^ = 0.1 
and Pi = 0.2. There is no peak splitting in Fig. (a) 
because most impurities are outside of the gray region 
in Fig. 13 On the other hand, the results in Fig. jSl (b) 
show the split of the ZBCP at T = 0, where we con- 
sider high density weak impurities with for ViNq — 0.02 
and Pi = 0.6. This is because most impurities are in- 
side of the gray region in Fig. [SI The splitting peaks 
merge into a single peak under finite temperatures such 
as T = O.lEzEP- In Fig.|Sl(c), we show the temperature 
dependence of the zero-bias conductance in Fig. (HI (b). 
The results show the reentrant behavior of the zero-bias 
conductance, which was found in the experiment. '^^ In 
Fig. El (d), the peak position {6eV) in (b) is plotted as 
a function temperatures. Since 6eV is about 0.15 Ezep 
at T = 0, peak splitting is washed out at high temper- 
atures such as T = 0.11 Ezep- High density impurities 
with weak random potential are responsible for the split 
of the ZBCP in low transparent junctions. 

Several experiments^Si^ show a sensitivity of the con- 
ductance peaks to external magnetic fields. Here we dis- 
cuss the conductance in the presence of magnetic fields. 
The effects of magnetic fields are taken into account phe- 
nomenologically by using the Aharonov-Bohm like phase 
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FIG. 6: The conductance in the presence of impurities dis- 
tributed randomly in the range of 1 < Xjkp < 20, where pi is 
the dimensionless area density of impurities near the interface. 
The conductance for low density strong impurities is shown 
in (a) with ViNo = 0.1 and pi = 0.2. The conductance for 
high density weak impurities with ViNo — 0.005 and pi = 0.6 
are shown for several choices of temperatures. The zero-bias 
conductance and the peak positions in (b) are plotted as a 
function of temperatures in (c) and (d), respectively. 



shifli^^^ of a quasiparticle. Since the impurity scatter- 
ing in magnetic fields itself is a difficult problem to solve 
analytically, we neglect the interplay between magnetic 
fields and impurity scatterings. Within the phenomeno- 
logical theory^ effects of magnetic fields is considered 
by replacing E in Eq. by £' + |Ao cos 7 sin 7 1 0b as 



(0) 

9 b 



7r/2 

d'y 

-7r/2 



Aq cos'^7 sin^7 
E]^z^ + Ag COS67 sin27 ' 



Eb =i? + 2Ao| C0S7sin7|(/)B, 
(f>B =271'— tan7 — Bq tan7, 



(56) 
(57) 
(58) 



where (po = 2'Khc/e and Bq = 1.0~^ corresponds to i? = 1 
Tesra. A quasiparticle acquires the Aharonov-Bohm like 
phase shift 0^ while moving near the NS interface. In 
a previous paper, we found that ZBCP in clean junctions 
remains a single peak even in the strong magnetic fields^"^ 
as shown in Fig.[7|(c), where zq = 10 and T = 0.05i?zEP- 
In Fig. (a) and (b), we show the conductance in the 
presence of low density strong impurities and high den- 
sity weak impurities, respectively, where Vi and pi are 
same as those in Fig.El(a) and (b), respectively. A tem- 
perature is fixed at T = 0.05£^zep- In contrast to clean 
junctions in Fig.[7|(c), the ZBCP in disordered junctions 
splits into two peaks under magnetic fields as shown in 
Fig. [3a) . The results within the phenomenological the- 
ory indicate that the sensitivity of the ZBCP to magnetic 
fields depends on the degree of impurity scatterings. In 



insets, peak positions {5eV) are plotted with circles as a 
function of magnetic fields. For high density weak impu- 
rities in Fig. 13(b), we also found that the degree of peak 
splitting increases with increasing magnetic fields. In the 
limit of the strong fields, 5eV tends to saturate as shown 
in the inset. These characteristic behavior are found in 
the experiment. "^^ Although we have well explained the 
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FIG. 7: The conductance under external magnetic fields for 
low density strong impurities with ViA'o = 0.1 and pi = 0.2 
are shown in (a), where T is fixed at 0.05i?zEP. Those for 
high density weak impurities with ViNo = 0.005 and pi = 
0.6 are shown in (b). In insets, peak positions are plotted 
as a function of magnetic fields. The conductance of clean 
junctions is shown in (c). 

characteristic behavior of the experimental conductance 
peaks under magnetic fields, the applicability of the phe- 
nomenological theory in the presence of impurities is still 
unclear. This issue would be addressed more clearly in 
an exact numerical simulation. 



IV. DISCUSSION 

In experiments, the split of the ZBCP has been re- 
ported in overdoped high-Tc superconductors. The 
heavy carrier doping may bring a number of defects or 
imperfections in superconductors. It is evident that the 
impurity scattering is unavoidable even in underdoped 
high-Tc superconductors. The split of the ZBCP would 
be found in underdoped superconductors if they have bad 
sample quality. In the same way, the split of the ZBCP 
would not be found even in overdoped superconductors 
if their sample quality are good enough. The sample 
quality is a key factor for the spilt of the ZBCP. This ar- 
gument is consistent with an experiment^ where the po- 
tential disorder is artificially introduced to the NS junc- 
tions by the ion-irradiation and washes out the ZBCP in 
the limit of strong disorder. 
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In order to compare the theoretical results in this pa- 
per with experiments, we may consider the impurity 
scattering in normal metals. The total resistance (R) 
in the dirty normal-metal / (110) d-wave junction can 
be described simply hy R ~ Rd + Rns!^^'''^ where 
Rd = ^/Gd is the resistance of the dirty normal metal 
and Gns = ^/Rns is the conductance discussed in this 
paper. The equation indicates the absence of the proxim- 
ity effect in the dirty normal metali2£*2& The height of the 
ZBCP can be reduced to reasonable values in the pres- 
ence of the impurity scattering in normal metals because 
the total conductance is given by 1/{Rd + Rns)- The 
degree of splitting depends on parameters such as the 
potential of impurities (Vj), the position of impurities 
(xj) and the transparency of the junction (Tb ^ l/^o)- 
In particular, Tb is the most important parameter. As 
shown in Figs. |31 and IHl the degree of splitting is roughly 
given by -Ezep = ^o/^o ~ ^qTb- Thus it is possible to 
choose Tb to fit the degree of splitting with that found 
in experiments. The amplitude of the pair potential is 
about 30-40 meV in typical high- Tc materials. The de- 
gree of splitting is then estimated to be 1.2-1.6 meV for 
zq = 5, which is almost consistent with that found in ex- 
periments, for instance, 2 meVi^ It is also necessary to 
consider electronic structures of high-Tc materials for the 
quantitative agreement of the splitting width in theories 
with that in experiments. 

In quasiclassical Green function theories, the conduc- 
tance is proportional to the density of states at the sur- 
face of superconductors. We find that E dependence of 
the density of states at the interface in Eq. is ap- 
parently different from that of the conductance even in 
the clean junctions shown in Eq. 1)46(1 . We also find that 
NsiE, X = 0) ^ No <S: NsiE,x = ^o) as shown in Fig. [3 
(b). The density of states near the interface averaged 
over ^0 in this paper may corresponds to the surface den- 
sity of states in the quasiclassical approximation, where 
the rapidly oscillating part of the wave function are ne- 
glected and the smallest length scales is given by ^o- It 
is impossible to directly compare the present theory with 
the quasiclassical Green function theories because the po- 
sition of impurities is a key parameter for the split of the 
ZBCP in our theory. 

The Abrikosov-Gor'kov (AG) thcory^^ is a useful ap- 
proach to discuss the impurity scattering in superconduc- 
tors. The applicability of the AG theory is limited to the 
superconductivity in diffusive metals, where the mean 
free path is much smaller than the size of the disordered 
region in superconductors. Here we briefly discuss a re- 
lation between the AG theory and ours. Since the AG 
theory assumes the diffusive transport regime, we define 
the mean free path of a quasiparticle in the Born approx- 
imation, 

i =VFT, (59) 

- =2^p,iVoF,2^ (60) 
r 

p. =N,/WL„ (61) 



where we consider W x disordered region on two- 
dimensional superconductor as in Fig. El The impurity 
density pi can be replaced by the dimensionless impurity 
density pi — {2Tr/kp)'^pi. The ratio of the mean free path 
and the coherence length is given by 



^0 ^^FP^{NoVi)^ ' ' 

In this paper, €/^o is larger than unity even in the limit of 
Pi = 1 because we assume Aq/pf = 0.1 and A^o^ < 0.1. 
On the other hand, the dirty limit is defined by a rela- 
tion £/£^o ^ 1. When we choose Ls being a few coherence 
length as shown in Fig. |^ where ~ S^o? we find that 
£/Ls is still larger than unity. The disordered region 
is not in the diffusive regime but in the quasi-ballistic 
regime because the diffusive regime is characterized by 
a relation £/Ls ^ 1. Thus it is basically impossible 
to apply the Abrikosov-Gor'kov theory to the model in 
this paper. The scattering theory^ used in this paper 
is a suitable analytic method to discuss the conductance 
of such NS junctions. Although the impurity scattering 
near the interface is very weak in normal states, it dras- 
tically affects the conductance below the critical temper- 
ature. Impurities located at the resonant states seriously 
suppress the degree of the resonance even if their poten- 
tials are weak. In the AG theory, the real part of the 
self-energy is usually neglected. In our theory, this ap- 
proximation corresponds to an equation S2 = 0. However 
■32 plays an important role in the peak splitting. When 
we omit S2, I2 and It, must be also zero, which leads 
to positive Tj and no splitting irrespective of ViNg and 
Xjkp. 

We show that the impurity scattering causes the split- 
ting of the ZBCP. The conclusion, however, does not deny 
a possibility of the BTRSS. In a previous paper )S we as- 
sume the s + id symmetry near the NS interface and nu- 
merically study the tunneling conductance. The results 
show that the split of the ZBCP is insensitive to the po- 
tential disorder. Thus peak splitting would be always 
expected in low transparent junctions^ if the BTRSS 
appears at the NS interface. At present, we have only 
limited information on the BTRSS within the mean-field 
theories. To understand the nature of the BTRSS fur- 
ther beyond the mean-field theory, we have to analyze 
electronic structure of high-Tc superconductors based on 
microscopic models and make clear effects of the surface, 
the electron correlation and the random potentials on 
the superconducting state. This is an important future 
problem. 

Since the formation of the ZES is a universal 
phenomenon in superconductors with unconventional 
pairing symmetries, the ZES is also expected at 
a surface of spin-triplet superconductors 1^ It is in- 
teresting to study effects of impurities on trans- 
port properties in spin-triplet superconductor junc- 

46.78.81.82.83.84.85.86.87.88.89.90.91.92.93.94 
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V. CONCLUSION 

We have discussed effects of impurity scatterings on 
the conductance in normal- mctal/rf wave superconduc- 
tor junctions. The conductance is calculated from the 
Andreev and the normal reflection coefflcients which are 
estimated by using the single-site approximation. We 
consider impurities near the junction interface on the su- 
perconductor side. The strength of the impurity scatter- 
ing strongly depends on the transparency of the junction, 
the position of impurities and the energy of a quasiparti- 
clc because the ZES are formed at the NS interface. We 
conclude that the impurity scattering causes the split of 
the zero-bias conductance peak. The results are consis- 
tent with previous numerical simulations. We have also 
shown that characteristic behaviors of the conductance 
spectra at finite temperatures and under external mag- 
netic fields qualitatively agree with those reported in ex- 
periments. 
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APPENDIX A: TRANSMISSION AND 
REFLECTION COEFFICIENTS 



In the clean NS junctions, the transmission and the 
reflection coefficients can be calculated from the appro- 
priate boundary condition of the wave function. The 
calculated results are shown below. 
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For instance, is the transmission coefficients from 

the electron branch in a superconductor to the hole 
branch in a normal metal. In above coefficients, we use 
a relation qf = kf' ~ ki for simplicity. The conclusions 
in this paper remain unchanged in this approximation. 
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APPENDIX B: GREEN FUNCTIONS 

The real space retarded Green function in clean junctions can be calculated by using the transmission and the 
reflection coefficients in Appendix A. For x < x' < 0, the Green function from a normal metal to a normal metal is 
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For X > x' > the Green function from a superconductor to a superconductor is 
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where if is the phase of a superconductor. For a; > > x', the Green function from a normal metal to a superconductor 
is 
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For X < < x' , the Green function from a superconductor to a normal metal is 
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